Microlocal analysis of scattering data for 
nested conormal potentials 



Suresh Eswarathasan 

Abstract: Working in the time domain, we show that the location of the singularities and 
the principal symbol of a potential that is conormal to nested submanifolds S'2 C 5*1 C , 
for n > 3, can be recovered from the backscattering as well as from the restriction of the 
far-field pattern to more general determined sets of scattering data. This extends the 
work of Greenleaf and Uhlmann where the potentials considered are conormal to a single 
submanifold S C M". We utilize the microlocal analysis of the wave operator D — df — 
and multiplication by a nested conormal distribution in order to study their action on 
spaces of conormal-like distributions. 

1. Introduction 

Consider tlie potential scattering problem for the wave equation: 

{dt - A + g)M = in R" X M, 

u = 6(t — X ■ oj) for t « —p, (1.1) 

where g is a compactly supported time- independent potential and u G S*""^ 
is varying. Here p is any value such that supp{q) C < p}. One obtains 
a scattering kernel for a fixed q, say in G C^(]R"), from the Friedlander 
radiation condition 11 81. 



n — 1 

aq{t,6,u)= \im r~dtu{t + r,r9, 00), (1.2) 

r— >-oo 



where (r, 6) are polar coordinates. The scattering map $ which sends q to 
aq is nonlinear and overdetermined and there has been much interest in the 
inverse problem of determining q from ag. Since ag is overdetermined, it is 
naturally also of interest to try to reconstruct q from the restriction of ag to 
various submanifolds of M x S*""^ x 5*""^. We write the scattering amplitude 
for the analogous stationary potential scattering problem as aq{X,6,u). 
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1.1. Statement of the problem and main result 

The class of g's to be considered in fll.ip are those that have singularities 



conormal [12| to a nested pair of sub manifolds 5*2 C 5*1 of M", denoted by 
(S'i,S'2), of arbitrary codimension. The inverse problem we solve consists 
of determining these submanifolds and the principal symbol of g, which is 
enough to determine the singularities of g, from the leading singularities of 
the backscattering a\^, where M = {9 = —uo} C M x S^~^ x 5""^, which is 
a distribution on M x S"""^. We treat similar determined sets of scattering 
data. 

It is shown that a, away from cj's that are tangent to either of the sub- 
manifolds, is the sum of a paired Lagrangian distribution associated to two 
cleanly intersecting reflected Lagrangians, two reflected Lagrangian distri- 
butions, and a single peak Lagrangian distribution, modulo Sobolev errors. 
Although the strongest singularity lies on the peak Lagrangian, as is well 
known in the physics literature, we show that it is the restriction of the re- 
flected Lagrangians and their points of intersection to various submanifolds 
of scattering data in M x S"""^ x S*""^ that determine the singularities of q. 
The precise theorem is the following: 

Theorem 1. Let 5*2 C 5*1 C M" he smooth nested submanifolds of codimen- 
sion di + d2 and di, respectively. Assume that q is compactly supported and 
is conormal to the nested pair (5*1, 6*2) of orders Mi and M2. Furthermore, 
suppose that 

M2 > -d2 and Mi < -rfi - y + 1 or 
M2 < -d2 and Mi < -di + 1, with 

+ ^ < mi{-^^^^{di + d2),-di-d2 + l}ifn>5, 
2 n 

and 

+ ^ < mi{-^^^^^,-di-d2 + l}ifn = 3andA. 

Then Si, S2, and the principal symbol of q are determined by the singularities 
of a restricted to the backscattering surface {6 = —u} C M x 5*"^^ x 5*""^. 

In fact, we prove a stronger result that shows that for a given sub manifold of 
scattering data with certain geometrical properties, the inverse problem can 
still be solved. See Section [H] for a more detailed statement. 
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The orders —^^(^1 + ^2) and in Theorem [T] are present so that we 

can follow the scattering theory of Lax and Phillips for short range potentials 
in certain L^{W^) spaces; see |26|]. The restriction of Mi + ^ < —di — d2 + l 
is needed in order to have a series describing singularities, which appears in 
a later section, stabilize in an appropriate sense. Moreover, the potentials 
under consideration are allowed to blow up. Specifically, for M2 < — ^2, q 
blows up on 5*1 and is continuous on S2, while for < —di — ^ + 1 ~ Mi < e 
and < M2 — d2 < e, q blows up on 5*2 and is continuous on Si\S2- Hence, 



there is no size restriction on q, in contrast to, e.g., [28 



Potentials that are conormal to a single submanifold were dealt with by 



Greenleaf and Uhlmann [10|; this paper closely follows the time-dependent 
approach taken in 10| and generalizes the results to nested q. Using Propo- 
sition 3.7 in ll| involving the intersection of classes of paired Lagrangians 



over all orders M2 , it follows that Theorem [T] covers the main result of [lO|. 

A significant difference between the work in this paper and that of [10|] 
is the new, more complicated geometry that arises when using an approxi- 
mation method, the understanding of multiplication by q on Sobolev spaces 
and other classes of distributions, and the appearance of distributions that 
are associated to cleanly intersecting triples and quadruples of Lagrangians. 

We conclude this section by noting that even for an arbitrary Lagrangian 
distribution u, calculating the blowup rates that assist in finding which 
space u belongs to is difficult without some additional assumptions on the 
Lagrangian. Hence, from the viewpoint of the Lax-Phillips scattering theory, 
assuming that u is in some conormal category is a reasonable restriction. 



1.2. Previous results 

1) Fixed angle scattering: 



One wants to determine information about the 



potential q from the scattering amplitude restricted to a fixed incident angle 



6*0 G S"""" . Stefanov [29[ proves uniqueness of the potentials under a small- 
ness assumption and Ruiz jisj shows that the Born approximation determines 
a "close" approximation of g G H^{MP) for n = 2 and 3. 
2) Fixed energy scattering: Here we set A = Aq in the scattering ampli- 



tude. Uniqueness and reconstruction results are obtained by Nachman [22 



Novikov 



24| ]. and Ramm 27| for dimensions n > 3. Sun and Uhlmann (32| 
show the recovery of singularities for n = 2. In this setting, the main tech- 



nique seems to rely on the complex exponential method used in [33| and the 



9-method utilized by Nachman 23 
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3) Backscattering: As mentioned earlier in Section [H we set 9 = —u in the 
scattering amplitude. Uniqueness under a smallness assumption is obtained 
by Lagergen 15|] for n = 3 and recovery of singularities for n = 2 by Ola et. 



al [25[. Generic uniqueness is proven in Eskin and Ralston 



1.3. Outline of the paper 

This work deals with the inverse scattering problem in the time domain 
and is done in order to utilize a more geometrical approach through the 
microlocal analysis of the wave equation and multiplication by conormal-like 
distributions, which allows precise tracking of singularities. 

Section [2]gives the necessary background from microlocal analysis with an 
emphasis on Lagrangian distributions, Fourier integral operators, and paired 
Lagrangians distributions. Here, we introduce the space of nested conormal 
distributions for which our potential q will be an element. 

In Section |3l we develop a new kind of Lagrangian-type distribution that 
will be seen in Sections El and [71 This work is a fairly straightforward 
generalization of the theory of nested conormal distributions. 

Sections HI contains a number of geometrical computations and lemmas 
necessary in solving our inverse problem. 

The main calculations in solving the forward problem begin in Sections [51 
and [HI where we try to understand the mapping properties of two operators, 
namely multiplication by q and the parametrix to the wave equation. Section 
[3 finishes the description of the leading singularities of the solution u to 
(11. ip using the Born series, relying heavily on the results from the previous 
sections. 

The inverse problem is solved in Section [HI where we describe the scatter- 



ing kernel ag using the Lax-Phillips scattering theory [16| and the Born series 
from Section [71 The paper concludes by showing that the various restrictions 
of the scattering kernel continue to determine the singularities of q. 

The paper is a revision of the author's Ph.D. thesis at the University of 
Rochester. 



2. Preliminaries 

Throughout this paper, unless otherwise specified, X will represent a 
smooth manifold of dimension n, {T* {X)\0 , ut* {x)) will be the cotangent 
bundle of X with 0-section deleted, equipped with the canonical symplectic 
2-form ut*{x) = ^id$,i A dxi. See for an overview of symplectic geometry. 
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The relation a < h denotes a < Ch for some constant C > 0, which may 
depend on some parameters but not those of interest. 

2.1. Fourier integral distributions 

Let S" C X be a smooth submanifold of codimension k. Suppose S = 
{x G X : hj{x) = 0, 1 < j < A;} is a local representation of S with {Vhj : 
1 ^ J ^ ^} linearly independent on S. 

Definition 1. We say that a distribution q is conormal to S of order /x, G 
R, zf 

q{x)= I e'^^^^^''^^^a{x-e)d9, (2.1) 

with a{x; 9) G ^^(X x (]R''\0)) and the space of all these is denoted by I^{S). 
We call a the symbol of q where \d2dga{x; 6)\ < {9)^^^"^ for {x,y) G K where 
the constant associated to < depends only on K,a, and 7. 

Here, we use the standard notation (9) = (1 + 16*^)^, the Japanese bracket 
of 0. 

If —k < ji < then g has a specific blowup rate: \q{x)\ < C-dist{x, S)~^~^] 



see Section 6, Section 4 of [SlJ for the calculation of this estimate. Examples 
of such g's are the surface measure dS along 5* as well as a function that has 
a Heaviside singularity across S. These distributions lie in I^{S) and I~^{S), 
respectively. 

For the conormal distributions q associated to a submanifold 5", WF{q) C 
N*{S) where N*{S) is the conormal bundle of the S. These, and many 
other kinds of geometric distributions, fall under the ubiquitous category 
of Lagrangian distributions; an important part of microlocal analysis is the 
study of such distributions. For the basic theory of conic and Lagrangian 
submanifolds of T*(X) and nondegenerate phase functions parametrizing 
Lagrangian submanifolds, see (l2| . 

Definition 2. Let A C T*(X)\0 be conic Lagrangian submanifold. The 
Hormander space /'"(A) of Lagrangian distributions associated to A consists 
of locally finite sums of distributions of the form 

u{x) = [ e'^("'^)a(x; 9)d9, (2.2) 



where (p is a nondegenerate phase function parametrizing A, a G S**""*"* 2, 
and WF{u) C A. 
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As mentioned previously, distributions conormal to a submanifold S fall 
into the class of Lagrangian distributions with A = N*{S). By definition, 



11^ (S) = /^-4 + 2(iV*(^)). The work of Melrose [19 gives an alternate char- 
acterization of the conormal distributions. 

Theorem. Let S G X be a smooth submanifold. The space of conormal dis- 
tributions on X with respect to S of order m, I"^{X] S), can be characterized 
as the set of all distributions u G T''(X) such that 

Vi...VkU e if~™-^'°°(X), 

where the Vj 's are C°° vector fields on X which are tangent to S , and 
H^'°°{X) denotes the Besov space of order s G M. 

We refer to this notion as the iterated regularity characterization of conor- 
mal distributions. A similar characterization exists for Lagrangian distribu- 



tions [13[. For more on Besov spaces, see [3C 

Now, let {T*{X)\0,cot*{x)), iT*{Y)\0,ijjT*{Y)) be the cotangent bundles 
of the smooth manifolds X and Y with the 0-sections deleted and respective 
symplectic forms. (T*{X)\0 x T*(Y)\0,ut*{x)xt*(y)) is a symplectic mani- 
fold with respect to the twisted 2- form ujt*(x)xt*{y) = t^l^t*(x) — 7'"R'^r*(y), 
where vr^ and tir are the left and right projections from our product space. 
If A c r(x)\o X T*(r)\o is a conic Lagrangian submanifold with respect 
to this 2-form, then A' = {((x; ^), (y; 77)) : ((x; ^), (y; — 77)) G A} is called a 
canonical relation ; A' is Lagrangian with respect to the twisted symplectic 
form t[Iujt*{x) - t^rUt*(y)- 

Definition 3. F G /'"(X, Y\ C) is a Fourier integral operator of order m 
(abbreviated by FIO) if the Schwartz kernel of F, KF{x,y), is an element of 
the space I"'{X,Y;C'). 

A natural question to consider: when is the composition of two FIOs 
again an FIO? 

Theorem 2. (Hormander 1971 ) Let Fi G I""' (X, F; Ci) and F2 G I""^ (Y, Z; C2) 
Suppose that x rtl T*{X) x AT.(y) x T*{Z), where rh denotes transverse 
intersection. Then Fi o F2 G /™i+™2(X, Z; Ci o Cg)- 

Operator compositions that fall under the hypothesis of this theorem are said 
to satisfy the transverse intersection calculus. 
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2.2. Paired lagrangians and nested conormals 

Classes of distributions associated with two cleanly intersecting Lagrangians 
manifolds were introduced by Melrose and Uhlmann and Guillemin and 



Uhlmann [llj in order to construct parametrices for systems of pseudodif- 



ferential operators that arose in various settings. For the purposes of this 



paper, we define the class through multiphase functions [20[ and symbol- 
valued symbols joj, as follows. 

Definition 4. J^] Let M and N he smooth submanifolds of X . Suppose that 
M n N is also smooth. Then for p ^ M (1 N , M and N are said to be 
cleanly intersecting at p if 

Tp{M) n Tp{N) = Tp{M n N). 

Moreover, {M, N) is a cleanly intersecting pair in codimension k if M and 
N are cleanly intersecting for all p & M (1 N and M (1 N is codimension k 
in both submanifolds. 

Definition 5. Let (Aq, Ai) be a cleanly intersecting pair of Lagrangians in 
codimension k in T*{X)\0. Suppose Aq G Aq n Ai andT C X x (]R^\0) x 
is an open conic set. A multiphase function cj) parametrizing the pair (Aq, Ai) 
is a function 9; a) G C°°{T) such that 

1. (f)o{x;6) := (j){x]6]0) is a nondegenerate phase function parametrizing 
Aq in a conic neighborhood of Aq, and 

2. (t>i{x; {9, a)) := 9; cr) is a nondegenerate phase function parametriz- 
ing Ai in a conic neighborhood of Xq. 



Example: The work of Guillemin and Uhlmann proves that any two 



pairs of cleanly intersecting Lagrangians are microlocally equivalent. One 
can thus consider the model pair (Ao,Ai) in T*(]R") where Aq and Ai are 
the conormal bundles of {x = (xi, ...,x„) = 0} and {x' = {xk+i, ■■■,Xn) = 0}, 
respectively, so that 

Ao = To*(R^) = {(0;0:eeM"\0} 

Ai = {{x, : Xk+i = ...=Xn = 0,^i = ...=^k = 0}. 

Then (p{x; 9', a) = x ■ {9', a), for {9', a) G (M"-^\0) x is an example of a 
multiphase parametrizing (Ao,Ai). 
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The singularities of distributions represented by an oscillatory integral is 
dependent on the interaction between its phase and symbol. As this discus- 
sion leads to a more general interpretation of Lagrangians distributions, it is 
natural to expect that the symbols themselves will generalize. 

Definition 6. The space S^'^'^''^{X x (M^\0) x M]"^) of symbol-valued symbols 
is the set of functions a{x;9;a) G C°°{X x M'^i x R'''^) such that, for every 
K ^ X, (a,/3,7) G Z^^ X x Z", the following differential estimate holds: 

\dld^,d'^a{x-Acr)\<{e,aY-\'^\aY-\P\ 
for X & K. If \6\ > \a\ on the support of a, we say that 6 is elliptic to a. 

Next, we define a generalized class of Fourier integral distributions associ- 
ated with a cleanly intersecting pair of Lagrangians. The following definition 



is a modification of the original formulation found in ^JO . 



Definition 7. Let (Aq, Ai) be a cleanly intersecting pair of Lagrangians in 
codimension k in T*{X)\0. The space of paired Lagrangian distributions of 
order p, I G M associated to (Aq, Ai), denoted by P'^Aq, Ai) , is the set of all 
locally finite sums of elements of P^^{Aq) + /*'(Ai) and distributions of the 
form 

u{x) = [ e*'^(^'^'")a(x; 9; a)d9da, (2.3) 



where a G S^'^X x (M^\0) x R''), with p = p + l + ^- I = -I - |, 

and (f){x] 9] a) is multiphase parametrizing (Ao,Ai) on a conic neighborhood 
of a point Aq G Aq fl Ai. 

Example: Consider u{x) = H{xi) ■ 5{x") where (xi,x") G M", H{xi) is the 
Heaviside function in xi, and 6{x") is the delta function. It follows that 
u G J?-i-t-5(Ar*({xi > 0}),N*{{x" = 0})). 

The potential in our scattering problem will be a similar kind of paired 
Lagrangian distribution. 

Definition 8. If S2 C Si G X are smooth submanifolds with codim{Si) = di 
and codim{S2) = di + d2, then N*{Si) and N*{S2) intersect cleanly in codi- 
mension d2- The space of distributions on X conormal to the pair {Si, S2), 
referred to as nested conormal distributions , is 

I^^^^^{Si,S2) = /"^^+^-^^+^-?'-*-^^HAr*(^i),iV*(^2)). (2.4) 
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If Si and S2 have defining functions {hi{x)}'f^^ O'^d {hi{x)}'f^'^^ , respectively, 
then u G /^^'^^(S'l, 5*2) can he locally expressed as 

[ e'P-^^(")^»]a(x;^',0")rfW, (2.5) 

where a e S^'^^'^'^^{X x (R"'i\0) xR'^^). 

Similar to the iterated regularity characterization of conormal distribu- 
tions, there exists one for the nested conormals. 

Definition 9. Suppose that Si and S2 are submanifolds of X that intersect 
cleanly. Then V(5'i, S2) is the space of vector fields that consists of those that 
are tangent to both Si and 5*2. 

Theorem. ^] The space of nested conormal distributions with respect to 
Si,S2, for 5*2 C ^i, denoted by I^^'^'^{Si, S2) , coincides with the set of all 
distributions u G such that for some Sq depending on Mi and M2, for 

all k>0, 

Vi...V,u G KUX) 

with Vj G V(S^ , 5*2) for I < j < k. Here, Hi^^iW") denotes the localized 
Sobolev space Jj/. 



3. Nested triple conormals 

In this section, we introduce a new class of distributions analogous to 
the nested pair conormals but with three Lagrangians intersecting pairwise 
cleanly with a smooth triple intersection. The work of Mendoza and Uhlmann 



21[ develops a symbol calculus for a class of distributions associated to an 
intersecting system of three Lagrangians. However the conditions imposed 
on their intersecting system are not satisfied in our situation, requiring us to 
take another approach. 

3.1. Iterated regularity characterization 

Definition 10. Suppose that Si, S2, and S3 are submanifolds of X that have 
smooth pairwise intersections, a smooth triple intersection, and are in gen- 
eral position; see f^, p. 83]. Then V(S'i, 6*2, 5*3) is the space of vector fields 
consisting of those tangent to Si, S2, and S^. 
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Definition 11. The space of triple product-type conormal distributions with 
respect to Si, S2, and S3, I{X; Si, S2, S^,), is the set of all distributions u G 
'D'(X) such that for some sq and all k > 0, 

Vi...V,u e H^^[iX) (3.1) 

with Vj G V(^i, ^2, ^3) for l<j<k. 

We will assume S3 C S2 C Si with codim(S'i) = di, codim(S'2) = di + ^2, 
and codim(S'3) = di + d2 + d^ . It is possible to introduce local coordinates 
in a neighborhood of S3 with 

51 = {xi = ... = XdA 

52 = {Xl = ... = X(ii+d2} 

53 = {xi = ... = Xd^+d^+ds}; (3.2) 
denote the points in M" as {x' , x" , x'" , x"") with x' = {xi, ...,Xdi), x" = 

-1-1, Xf^ -|-c(2 ) ! X " = {Xd,+d2+l, ■■■,Xd^+d2+d:i), andx"" = {Xd,+d2+d3+l, ■■■,Xn). 

Let (^', ^"") be the dual variables in phase space. The following is 



based on the iterated regularity characterization for nested pairs in [19 
also [9J. 

Proposition 3. If u e V'{W) and S3 C S2 C Si C R" are as m / f5^) . then 
u G /(M"; Si, S2, S3) if and only if there exists an sq eM such that 

D:,Df„D:,„D2""((^')''(^")'(^'")^«) e //^«(M"), (3.3) 

for all multi-indices a, (3, rj, 7, p, 6, a where \p\ > \a\, \6\ + \p\ > \a\-\-\f3\, \6\ + 
\p\ + \a\ > \a\ + \(3\ + \ri\. 

Proof : We find a local set of generators, over C°°(M"'), for the algebra 
V{Si, 5*2, 5*3) consisting of the differential operators that appear in (13. ip : we 
will do induction on the order of the operator in (13. 3p . 

Take V = ^tMx)D^. + ^fJd%M^)D., + ^t'd%'t%Mx)D^k 
+'^e=di+d2+da+iMx)Dx^, with ai,bj,Ck,de G C°°(M"), then V being tangent 
to 5*1, 5*2, and S3 means that the = at x' = 0, = bj = at x' = x" = 0, 
and Qi = bj = Ck = at x' = x" = x'" = 0, respectively. Then, 

XiD^^ for i < di,j < di, 

XiD^- for i <di + ^2, di < j <di + ^2, 

XiDxj for i < di + d2 + (is, di + d2 < j < di + d2 + c?3, and 

D^^ for k> di + d2 + ^3, (3.4) 
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gives our desired basis. 

Assume fl3.3p holds for |a| + + + I7I < P, where the induction step 
will be done on p. The application of an operator Vp+i G V(S'i, 6*2, 5*3), which 
is a linear combination of the vector fields in (13 ■4p . to the expression in fl3.3p 
gives 



= ^<^'X^D:,dI,,D:,,,DI^^^^^^^^ (3.5) 

+ ^h':;;:;mD^^,,D:.^^^^^^^^ (3.6) 



where j, k have the restrictions indicated in (13. 4 p and the coefficients are in 
C°°. Further commutation of the differential operators with the coefficients 
x', x", x'", and 0j gives the expression in (13. 3p with |a| + |/3| + |'7| + |7| < 
p + 1. ■ 



Proposition 4. // n G /(X; ^i, ^2, Sg), i/ien lyF(M) C Ar*(^i) U Ar*(52) U 
N*{S^). 



Proof : Using Definition 13.11 and the vector fields computed in the proof of 
Proposition [3l it follows that u is not singular in the x"" variable. ■ 

The iterated regularity definition is not limited to the configuration of 
nested submanifolds. The computation of a basis for V(S'i, 5*2, 5*3) for a more 
general triple needs additional assumptions and is more complicated. This 
is left to the interested reader. 

3.2. Oscillatory representation 



Definition 12. The space 5*^i'^^2,Af3(Xx (M^\0) xM^^ xM'^-') of triple symhol- 
valued symbols is the set of functions a{x; 6; a; r) G C°°(X x M^^ x M^^ x M^^) 
such that, for every K ^X, (a, /3, r/, 7) G Z^^ x Z^^ x Z^^ x Z\, the following 
differential estimate holds: 



for X G K. 

Proposition 5. Let u G /(M"; Si, S2, S3) with the local coordinates in ^3.2i) . 



Then 




(3.7) 



where a G 5'A^i.*^2,M3(]Rn x (Mf\0) x M'^^ x R'^s). 



11 



Proof : Following Melrose's original line of argument, also found in [qI, we 
can assume that u is compactly supported; otherwise, we can multiply by a 
smooth cutoff supported near 5*3. Take the Fourier transform in {x',x",x"') 
and using fl3.3p we get 

{e,e,CrdededCdx"") (3.8) 

for some sq when \p\ > \a\, \6\ + \p\ > \a\ + \f3\, \6\ + \p\ + \a\ > \a\ + \f3\ + \r]\. 
The Sobolev embedding theorem |7j tells us a(x; ^"') := u{^', x"") 
satisfies a triple symbol-valued symbol estimate for some Mi,M2, and M3 
depending on the dimension of our submanifolds Si, S2, and S3. ■ 

We will now generalize the proof of Proposition 1.20 in joj to show that 
distributions of the form (13. 7p are independent of the choice of coordinates. 

Suppose u is of the form (13. 7p with an integral that is absolutely conver- 
gent; otherwise, integrate by parts to lower the order of the triple symbol- 
valued symbol. A change of variables that preserves (13. 2 p must be of the 
form 

= Ejii Aij{y), l<i<di 
x^ = Tf■=l"B^^{y), di<i<di + d2 (3.9) 



= E ■it"'"""' C',,(2/), di + d2<l<di+d2 + d: 
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with Aij,Bij, and Cjj being smooth. Plugging this new coordinate transfor- 
mation into (13. 7p and setting 



= Ej=l ^ji^j + J2j=di+1 Bji^j + J2j=di+d2+l ^Ji^J^ 1 < ^ < C?l 

= Ejldi+i Bji^j + Ylt=dt+dt+i di <l <di + d2 

'^i = T.jLd2+d3+l ^ji^V di+d2<i <di 

+d2 + d^, 

(3.10) 

gives the new oscillatory representation of 

u{y) = J e^(s'-y'+^"-y"+^"'-y"%{y; E'; E"; E"')dE' dE" dE'" . (3.11) 

Here, the interested reader can prove h = a{x{y),C,'{y, E', E", E'"),C^"{y, E', E" 
,E"'),^"'{y,E',E",E"')) X |-^| is another triple symbol-valued symbol. This 
shows that (13. 7p is independent of the choice of coordinates. We can now 
make the following definition. 
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Definition 13. Let S3 C S2 C Si be submanifolds of X with codimensions 
di + d2 + 4, di + c/2, and di, respectively. Then lMuM2Mi(j^n. 
the space of locally finite sums of distributions of the form 



where a e S^^^^^^^^{R''x{Rf\0)xR'^^ xR'^'^) and {Hi{x)} , {Hi{x) = H2{x)}, 
and {Hi{x) = H2{x) = H^^x)} are defining functions for Si,S2, and S3, 
respectively. 



Proposition 6. Letu G S-^^ S2, S3). Then, away from N*{Si)n 

N*iS2)nN*iS3), 



for a e ^A^fi.^^2.A^3_ Therefore, 

N*{Si) = {(0,x",x'",x"";r,0,0,0):r^O}, 

N*{S2) = {(0,0,x'",x"";e',r,0,0) : {^,0 + 0}, and 

N\S3) = {(o,o,o,x"";e',e",r,o) : (e',r,n ^ o} (3.15) 



N*{Si) n N*{S2) n N*{S3) = {(0, O, O, x""; e', 0, 0, 0) : 0}. (3.16) 

Away from (KWl . we have {x',x",x"') ^ 0, C 7^ 0, or C" ^ 0. In the 
first case, integration by parts in the or C,'" variables contributes a C°° 

function, an element of I^^^{Si), or an element in /^^^'^^^(S'l, 6*2). 

Now, suppose C 7^ 0. Then either {|^"| > |^'|} or | > |^'"|}. Assum- 
ing the former case, we get that u e lMi+M2+hh+-'^+Y'''' -*^3-^ (A^*(5'2), 
N*{S3)) since a G 5'^^i+^^2,Af3_ j^y^Q being away from the intersection de- 
scribed in (13. 161) . ^' is allowed to be 0. As ^' cannot be elliptic to all the other 




(3.12) 



Proof : By using Definition [131 we can assume 




(3.14) 



which shows 
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variables, for otherwise we would be localized to the intersection in fl3.16p . 
it follows that we are either again in the region {|^"| > |^'|} or in the region 
^ l^'l}; the second region is discussed next. 
Consider C" ^ 0. Then either > or > The former 

case again shows that 

while the latter case implies 

as a G 5'^^i'^2+Af3_ Eg^Q]-^ Qf paired Lagrangian distributions found falls 
into a class described in the statement of the proposition. ■ 

4. Lagrangian submanifolds and canonical relations 

Each of the subsections below will be focused on the specific geometry in- 
duced by a given operator. In each case, we give the corresponding canonical 
relations for the operator, analyze the compositions with certain Lagrangian 
manifolds, and state properties of the resulting manifolds. 

4.I. Geometry determined by a multiplication operator 
Recall for q G C^(M"), multiplication by q, which we denote by Mq, is a 
pseudo differential operator of order 0. Moreover, Mg for q conormal of order 
H is an example of, in the language of [9l| , a pseudodifferential operator with 
a singular symbol of order and fi when using conormal notation. This 
notation will be described further in a later section. 

For the operator Mg with q G J*^^'*^^ (Si, 5*2), there are three canonical 
relations in T*(M" x R x 5""^) x T*(R" x M x S"-^) associated to it: 

Ar*(R"xRx5"-i)xT*(R"xRxS'"-i)) (4-1) 

Cs, = {{{x,t,u;;^,r,n),{x,t,u;^ + D,{{hi)ti)-9,r,n)) : (4.2) 
{x,t,uj;^,T,n) G T*(R" X R X 5"~^)|5,, and 6 G M.'^'\0}, 

where {/ij}*^^ defines Si and 

Cs, = {{ix,t,oj;^,T,n),{x,t,oj;^ + D,i{h,)tYn-0,T,n)) : (4.3) 
{x,t,u;^,T,Q) G T*(M" X M X S"-^)|52, and 6 G M'^^+^^VO}, 
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where {/iij^lf^^ defines 82- For our purposes, the three most important La- 
grangians in T*(R" x R x S"'^^) that interact with the canonical relations of 
Mq are 

A+ = {(x, X ■ — acj, (T, — cri* (x)) : (4.4) 

(x, t, w) G X M X a G M\0} 

Ai = {{y,t,oo;iy,0,0):{y,u)eN*{Si)} (4.5) 
A2 = {iy,t,u;u,0,0):{y,u)eN*iS2)}. (4.6) 

It is important to point of out these three Lagrangians are also the conormal 
bundles of 5+, ^i, and S2, respectively, where 5*+ = {t — x-w = : (x, t, w) G 
M" X M X 5"-!}. As Si ^ 5+ for i = 1, 2, it follows n 5+ = Si+ are smooth 
sub manifolds of M" x M x S"^^. Composition of the canonical relations in 
M and (gSI) with A+ give 

Ai+ = {{y,y ■ uj,u;iy - au,a,-(Til{x)) : (4.7) 
(y, u) G A^*(5i), a; G a G M\0}, and 

A2+ = {{y,y ■ u,u;i^ - auj,a,-(ril{x)) : (4.8) 
(y, z/) G Ar*(52), u; G a G M\0}. 

Once again the Lagrangians submanifolds described in 04.71) and (14. 8 p are the 
conormal bundles of 5"!+ and 6*2+. From this fact, it follows immediately that 
(Ai_|_, A2+) is a cleanly intersecting pair in codimension d2- Similar reasoning 
can be applied to other pairs of the above Lagrangians. 

4-2. Geometry determined by the wave equation 

The operator D^^, to be elaborated upon in Section 7, is another example 
pseudodifferential operator with a singular symbol. More, specifically D"^ G 

3 1 

I 2' 2 (Aj'.(KnxMx5"-i)xT*(K"xMxS'"-i)) C'n)) where 

Co = {(^ix,t,u;^,\^\,n),ix+it-s)-^^,s,u;^,\^\,n) 

(x, t,u) eW xRx s G M - 0, ^ G M", 

n G T:(5"-1)} (4.9) 

is the fiowout relation generated by the characteristic variety 

char{D) = {(x, t, co; ^, r, fi) : |r|2 = (4.IO) 
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The application of this canonical relation to certain Lagrangians yields 
new Lagrangians away from a prescribed set. As Aj, for i = 1,2, does not 
meet char{n), it follows that Cn o Aj = and any application of D"^ to 
distributions with WF not supported there will give a C°° function. Also, 
since A+ C char{n\), Cu o A+ = A+. 

Let 

= Ai+n char{D) 

= {|z/i|2-2z/i-a; = 0:(yi,z/i)GA^*(5i),a;G5"-i} (4.11) 

and 

= A2+nc/iar(n) 

= {W2\^-2u2-uj = 0:{y2,i^2)eN*{S2),ueS''-^}. (4.12) 

It is from this set that our new smooth Lagrangians will be flowed out. 
However, there exist points for which these fiowouts are not guaranteed to 
be smooth, therefore demanding that we impose some further restrictions. 

From (14. lip and (I4.12p . if z/i ■ = or z/2 ■ = 0, then Ui or 1/2 is forced to 
be therefore putting 0-sections into Ai or A2, which violates the 0-section 
assumption. Notice that this dot product is when u is tangent to the 
submanifolds 5"!+ or 5*2+ (In [lo|, such elements of S"""^ are referred to as 
tangential rays). Hence, it is natural to consider the open region where u is 
not tangent to 5*1 or S2. We define 

El := E^ n {z/i ■ a; ^ : {y,, u,) G N*{S,),u G 5"-^} 

and 

E2 := E^ n {z/2 • u; ^ : (7/2, 1^2) G N*{S2),co G 5""^}. 
After solving our characteristic equation, we arrive at 

51 = {(yi,yi ■^,w;z/i -cra;,cr, -cri* (a;)) : (2/1,^^1) e A^*(5'i), (4.13) 

\j. 12 

CO G = ^^}, and 

52 = {{yi.Vi- uj,uj]U2- (TUJ,(T,-ail{x)) : {y2,V2) & N*{S2), (4.14) 

u e b ,cr = 1. 

2z/2 ■ u 
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This reparametrization is done to show the above sets are smooth submani- 
folds of Ai+ and A2+, respectively. 

The Hamiltonian vector field associated to the symbol of □ is 

H --e — + r— 

dx dt^ 

where x and t are the spatial and time coordinates on T*{W^ x M x S"""^). 
In order to guarantee our fiowouts are smooth Lagrangians, we must start at 
points of El or S2 at which Hjj is transverse. In order for this to happen, 
cannot be in the tangent space of Sj, for i = 1,2. For if it were, then — ^ • ^ 
will lie in the horizontal tangent space of the manifolds parametrized in f l4.13p 
and fl4.14p . But under these parametrizations, —C, = —{v — au) and t = cr. 
The only way this can occur, in Ei for example, is if = and u lies in the 
tangent space of Si+, with either violating the conditions in f l4.13p . Hence 
Ha iti Sj for i = 1,2 and generates the following Lagrangian manifolds: 

Ai+ = {{y-r{iy-au),y-UJ + r,u;iy-cxuj,a,-aC{y)) : (4.15) 
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{y, v) e N\S{),uj e S"~\ r e M - 0, a = ^f^} 

A2+ = {{y-A^-^^)^y-^^'r.^\'i^-o^.o.-^^*u,{y))'- (4.i6) 



\v 
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(y, v) G iV*(52), a; G S^-\ r G M - 0, a = 



As all of our parametrizations come from a single coordinate patch, we 
can set coordinates equal to one another, solve the resulting equations, and 
find parametrizations for the intersections. Doing this for A+ and Ai+ shows 

u = 0. Similarly, when looking at the parametrization of A^_,_, evaluating the 

1^12 

expression at z/ = makes cr = 0, violating the 0-section assumption 
once more. This shows A+ fl Ai+ C The exact same reasoning gives 

A+nAi+ c s2\S2. 

The differential of the projection vr from the A^_,_ onto the spatial variables 

i*{In~di+ rdy{u - au)) v - auj rK*de{y - auj) ril{d^{auj))\ 

fiu-rdy{a)) -a -rK*de{cj) il{y) (4.17) 

h-i j 
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Figure 14.21 A schematic representation of the Lagrangians; the cone 

represents char{D). 



where j : Ty{Si) ^ Ty{W), i : T^(5"-i) ^ T^(M"), and k : c/i-plane ^ 
Ty(R"') are the differentials of the inclusions of the respective submanifolds. 
By using Euler's identity for homogeneous functions, it follows ^ G span 



rank of this projection is locally constant and equal to 2n — 1 if the rank of 
the submatrix 



is locally constant and equal to n — d. This will hold away from the lower- 
dimensional set {(— r)"~'^Met(— (iy(z/ — cru) — \ln-di) = 0}, which is the 
singular set for the top submatrix in (14.181) . We remind ourselves that r 7^ 
for otherwise our parametrization drops rank. Notice this determinant gives 
a polynomial in r, having only finitely many solutions in this parameter. The 
same arguments are used to show that the rank of vr from A.2+ is also 2n — l. 
We use Prop. 3.7.2] on conormal bundles to help conclude 

Lemma 7. Away from a lower-dimensional set in phase space, A"_^ and Ag.,. 
are the conormal bundles of hypersurfaces m M" x M x S^~^ . 

This fact will play a substantial role in the approach we take when solving 
our desired inverse problem, particularly in Section 8. 





(4.18) 
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We will now recall a simple fact from linear algebra. Suppose that Wi 
and W2 are subspaces of a finite dimensional vector space V. We know 
dim(Tyi + W2) = dim(iyi) + dim(iy2) - dim(iyi fl W2). If we consider two 
submanifolds M and iV of X, and assume p & M (1 N with the intersection 
being smooth. As Tp{M n N) C Tp{M) n Tp{N), it follows 

dim(TpM + TpN) < dim(TpM) + dim(TpA^) - dim(Tp(M n A^)). 

For the given flowouts, an upcoming calculation shows A"_,_ fl has dimen- 
sion 2n - d2. For A G A^^ n A^_j_, dim(TAA5'_^ + TaA^^) cannot be greater 
than 2n + (^2. Hence, if we can prove this sum of tangent spaces has dimen- 
sion at least 2n + c?2, then dim(Tp(M n A^)) = dim(Tp(M) n Tp{N)) and the 
intersection is clean. 

Let us label each parameter that appears in fl4.16p with the symbol 
above it. Also, let 6' G M'^^yO parametrize the fibers of A^_^ and {9' ,6") G 
]^di+d2^Q parametrize those in A2+. Observing the spatial spherical and r co- 
ordinates imply that a = a and u = the equality in the ^-coordinates tells 
us that V = V and that the fibers which are characteristic in A^*(S'i+) must 
coincide with those that are characteristic in A^*(S'2+). Furthermore, after 
setting the x coordinates equal and using the facts above, y = y. It is worth 
noting a = a is actually an implication of N*{Si+) intersecting A^*(S'2+). 
In conclusion, it follows that in either parametrization of our flowouts, the 
intersection is obtained by simply restricting ourselves to 5*2+ in the (x, t) 
coordinate or to the fibers of N*{Si^) in the (^, r) coordinates in the case of 
A°_,_ or Ag.,., respectively. As this restriction only drops rank by d2, A^_^nA2_,_ 
has codimension d2 in both flowouts. 

To show the dimension of the tangents do add up correctly, we con- 
catenate the differentials of the parametrizations in fl4.15p and (14.161) in the 
following manner: 



where the left submatrix corresponds to (14.151) and the right submatrix cor- 
responds to (I4.16p . We will first look at the concatenation of two submatrices 
in the above matrix, namely ^Xfn^^ and , which is 



f n*dei{u — au) K*d0,{v — auj) K*d0i,{v — auj)\ (aocW 
\ K*de'{-cruj) K*dg,{-aoj) K*dg„{-aid) J \ ■ ) 

Doing a row reduction by subtracting a multiple of the r-row from the ,^-rows 




(4.19) 
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gives us 

/ K*de'iiy) K*dg,{iy) K*dg.{u) \ 

\K*de'{-au) K*dg,{-au)) K*dg„{-aoj) J \ ■ ) 

As z/ = z/ on the intersection of our flowouts, the first di columns are 
equal to the second di columns of this matrix. Moreover, as u forms a frame 
in the fibers of N*{S2+), the last d2 columns are independent from the first 
di. Hence, the rank of this concatenation of submatrices is di + d2- The 
rank of the differential in (14.151) is invariant meaning that if we do the same 
row reductions on the submatrix -%||7y^, then the rest of the rows in the 

differential ^^^^y'^g^'^'^^ will always give us rank 2n — di. Now, for every 
rank of di given by this first submatrix, we get another d2 for free by our 
observation. Therefore (14.191) has rank at least 2n + d2 and A°_,_ fl A2_^_ is a 
clean intersection. 

It is crucial for the intersection of the flowouts to be nonempty even 
after we avoid our bad set S^\S2 U S^\S2 for otherwise the standard mi- 
crolocal analysis cannot be applied. This can be checked by examining the 
parametrizations. 

Remark 1. The above calculations imply the intersection of the flowouts 
determine one of the initial surfaces, namely 5*2. 

Remark 2. While A^_,_ and Ag.,. are vector bundles over hypersurfaces, the 
clean intersection of these Lagrangians implies the underlying surfaces inter- 
sect tangentially. 

4-3. Geometry associated to Radon transforms 

Let R : ^'(R") £'{R x ^""i) be the Radon transform 

{Rf){s,e,t,uj)= [ f{x,t,uj)da{x) (4.22) 
Jxe=s 

where da is normalized Lebesgue measure on the hyperplane {x-6 = s}. Ris 
an elliptic Fourier integral operator [3], with R G I~^{Cji) where Cr is the 
local canonical graph inside of T*(R x ^""^ x M x 5"^^) x T*(R" x M x S""-^) 
given by 

Cr = {[x ■ 9,9,t,u;a,ail{x),T,Q), (^x,t,Lo,o-9,T,Q) : 

{x, t, u) e R"+^ X e e (T G R - 0, r G R, 

Q G T*{S''-^)}. (4.23) 
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An operator F whose canonical relation is Cn will be used when construct- 
ing the scattering kernel in Section [SI When we apply Cr to A_|_, A"_^, and 
Af.,., remembering that the application of canonical graphs to Lagrangians 
always satisfy the transverse intersection calculus, we get 

C^oA+ = {{yuj,uj,yuj,uj;a,-ail{y),-a,ailiy)):yeW,{4:.24:) 

U {{-y- CO, -u, y-u,u; a, -(riUv), cr, -o"C(^/)) : V ^ I^". 
CO G G M-0}, 

CnoAU = {(r-,.i^^^,-(^^^,,.a; + r,a;;a, (4.25) 
ai\_^^{y - r{v - auj)),a,-ail{y)) : 
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{y,v) G N*{Si),uje S''-\r G M, a = with z/ ■ w ^ 0} 

21/ ■ 00 

(I (v — OIjO) 

U |(— T + ?/ ■ , z/ — crw, y ■ + r, (T, 

cr 

(yj* j.-au.) {y - r{u - ouj)), cr, -ailjy)) : 

(l/,z/) G N*{Si),coe S"-\r G R,a = with z/ ■ w ^ 0}, 

• 00 

C.oM, = {(r-,.<^i^,-<^i^,y.<. + r,<.;a, (4.26) 

gj* (.-.c.) - r(z/ - (To;)), g, -ffz* (y)) : 
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{y,v) G A^*(52),u; G ^"-\r G M, a = -^^^^ with z/ ■ w ^ 0} 

2z/ ■ 00 

(u — au) 

U |(— T + ?/ ■ , u — (Jco,y ■ CO + r,co; a, 

cr 

ai*_f^^{y - r{u - auj)), a, -ai*^{y)) : 
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(l/,z/) G iV*(52),w G S"-\r G M,(T = with z/ • w 0}. 

2z/ • CO 

We note that as Cr is a local canonical graph, i.e. locally the graph of a 
symplectomorphism, {CroAI^, C/joAg^) is a pair of Lagrangians intersecting 
cleanly in codimension 
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4-4- Geometry associated to pullbacks 

For to » 0, the mapping : M x 5""^ x S""'^ ^ M x S""'^ xRx ^""^ 
given by pt^ = (s, 9, u) = (to + s, 9, to, w), induces the restriction mapping 

p*, : V'p{R X 5""^ X M X 5"-^) ^ r''(M x 5""^ x 5""^), (4.27) 

where is the space of distributions whose wavefront is disjoint from the 
normals of ptg. We use this space in order to make this restriction map well- 
defined. We drop the subscript to from here on, as it is understood what pto 
means. It follows that p* is Fourier integral operator and p* G I^{Cp), where 

Cp = {(s, 0, w; r, $, fi), (to + s, 0, to, r, $,?7, i7) : s G M, and w G S""-^, 
(to + s, 0, to, u; r, $, r/, fi) G T*(R x S"^^ x R x S"-i)\0} (4.28) 

It was proven in that the compositions of Cp with the Lagrangians in 



fl4.24p . 04.251) . and 04.261) are all transversal. In fact, the second components 
in 04 241) . 047251) . and fl4726D vanish after the application of Cp. 

Definition 14. The peak Lagrangian is defined as 

A+ = CpoCRoA^ (4.29) 
= {iy,u,u;a, -<(!/),<(?/)) : y G M",u; G G M\0}. 

Definition 15. The reflected Lagrangians are defined as 

AU = CpoCnoAt^ = {{-y.{^-^^^+u),-^-^^^ 
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and 



r G M, 0" = with z/ ■ w 7^ 0} 

Iv ■ UJ 



= C,oC^,oA^, = {(-,.(i^^^+c.),-l^^,c.;(4.31) 

-azt._..)(y),<(y)) : {y,u) G N*{S2),u e S^-\ 



|z/R 

r GM^a = with u ■ u ^ 0} 

2v ■ tu 
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A direct calculation shows A"^nA2_,_ is a smooth submanifold of codimension 
d2 inside each of the reflected Lagrangians. 

The proof that A°^_ fl A2_^_ is clean is almost identical to that which shows 
A"_,_ n Ag.,. is clean. The concatenation of the differentials of the parametriza- 
tions in fl4.30p and f l4.3ip will have rank 2n + d2 if the concatenation of the 
submatrices ^^^gi^^ and j^'^gf'g - ,}-^ has rank di + d2- This follows immediately 

after realizing that '^^j^;^^ is just the composition of from (14. 7p . which 
has rank di, with the polar coordinate map, which has an invertiable differ- 
ential away from the origin. The same holds true for ^^^fp-^ in relation to 

^^^g/g),^ from (14.81) . which has rank di + d2. Moreover, has as its columns 

the first di columns of j^r^- The remaining reasoning is similar. 

Remark 3. Once again, the clean intersection calculations indicate that the 
reflected Lagrangians should give us substantial information about Si and 82- 
This will be justified in the final section of this paper. 

5. Multiplication by a nested conormal 

The operator Mq, which is multiplication by g G /^^'^^(S*!, 5*2), has wave- 
front set in the three canonical relations 

^r*(R"xMxS"-i)xT*(R"xRxS"-i)) (5-1) 



Cs, = {{x,t,u-^,T,n),{x,t,co-^ + D,{{h,)ti)-9,r,Q): (5.2) 
{x,t,u;^,T,n) e T*{W X M X 5''~^)|s,, and 9 e M"'^\0}, 

and 

Cs, = {{x,t,u;;^,T,n),{x,t,co;^ + D,{{h,)tV')-9,r,n): (5.3) 
{x,t,uj;^,T,n) E r*(M" X M X S''-^)\s„ and 6 E W^'+'^'\0}. 

An operator of this type no longer falls into the class of operators whose 
Schwartz kernels are paired Lagrangians; in fact, the above canonical rela- 
tions form a triple of pairwise cleanly intersecting Lagrangians. The lack 
of a developed theory for the compositions of such operators requires us to 
take some different approaches when analyzing its mapping properties on 
distributions. 
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Mg is technically defined on the space of distributions V that has the 
property {Cs^ o WF{v)) n = and {Cs^ o WF{v)) n = 0, for all v eV. 
However, the multiplicative results and Sobolev mapping properties of this 
section circumvent this technicality. 

5.1. Action on spaces of distributions 



The following is a result from [10 



Lemma 8. Let X he a manifold and Y.,Z <Z X he suhmanifolds with Y &i Z, 
then 

/^(r)-F'(Z) C I^''''' {¥,¥ n Z) ® {Z,Y n Z). (5.4) 
Moreover, if u E I^{Y) satisfies supp u gY, then 

u- I^{Z) C I^'^'{Y,YnZ). (5.5) 

We need another lemma stating the multiplication properties of distributions 
associated to nested pairs of submanifolds. 

Lemma 9. Suppose Y^, Yi, and Y2 are suhmanifolds of X such that Y2 G Yi, 
F+ &\ Yi, and F+ ^ K,- U Ui GG !"{¥+), U2 G /^^^■^^^(Fi, F2) with supp Ui 
C y+ and U2 supported microlocally near N*{Yi) fl N*{Y2), then 

U1U2 G n r+, Y2 n r+). (5.6) 

Note that this result is a direct analog of the previous lemma for nested 
conormal distributions. 

Proof : We can introduce local coordinates (x', x" , x'") G M!^^ x'^'^Wx M"-'^i-'^2 
near the point Xq G Y^nY2 such that Xq = 0, F+ = {x' = 0}, Yi = {x" = 0}, 
and Y2 = {x" = x'" = 0}. Under these assumptions, 



u^(x) = / e'^^'<"^a{x;e)d^' 
for a G S^{X x (R^^^yo)), and U2 as the representation 

U2ix)= [ e^(^"-«"+^"'-«"')6(x;r,r)«'" 

for b G 5*^1 ■^^^(X X (]R^i\0) x R'^^)^ Multiplying the distributions leads us to 
the oscillatory integral 

{u,u2){x)= I e^(^''«'+^''-€''+^'''-«''\(x;eoK^;r,r)^e'«'''. (5.7) 
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since U2 being supported near the aforementioned intersection is equivalent 
to having b locahzed to the region = {^"| > |^'"|}. There are three regions 
in which can he: I = > > ICl}, H = {1^1 > W\ > 1^1}, and III 
= {W\ ^ W'\ ^ Computing the size of the Japanese brackets along 

with their orders, it follows that regions II and III correspond to the nested 
conormal I{Yi, Y2nY+) and the nested triple conormal I{Yi, YiHY^, F2nF+), 
respectively and of certain orders. A general fact involving the multiplication 
of distributions from [12[ says that 

WF{uiU2) C N*{Y+) U N*{Y+ f] Y^) U N*{Y+ n Y2) 

as supp Ui C y+. All three terms above lie inside the class J/^'^^i'^a^y^^ Yi n 
F+, Y2nY^), defined in Section [31 Moreover, the region in which this integral 
has the singularity corresponding to the triple intersection is I because of ab 
satisfying a triple symbol- valued symbol estimate, namely ab G S^'^^^'^^'^. ■ 

5.2. Sobolev estimates 

In this section, we prove estimates for the operator that multiplies by 
a nested conormal distribution. Previous estimates for multiplication by 
a standard conormal distribution were obtained in jof to describe operators 
with two canonical operators, one being a fiowout relation and the other being 
the diagonal relation. This work relied heavily on a composition calculus for 
certain kinds of paired Lagrangians developed by Antoniano and Uhlmann 
[l[ and follows Hormander's method for obtaining estimates for Fourier 
integral operators associated to canonical graphs. 

As stated in Section HJ Mq is now associated to three canonical rela- 
tions, two of which are fiowout relations. The lack of a formal composition 
calculus for operators associated to three intersecting Lagrangians forces us 
to take another approach. We will make use of an observation of Melrose 
(lof (see also joj) that elements of J^^^'^^^j^y^ y^) can be decomposed into 
a sum of two conormal distributions with (|, |) symbols; this is called a 
parabolic decomposition. From here, we use an orthogonality argument to 
essentially retrieve the same estimate as for multiplication by a standard 
conormal distribution, but in the (|, |) case. For more on (|, |) symbols, see 



3ll . Chapter 7]. 



Proposition 10. Let Yi C Y2 C X with dim X = n, codim Yi = di and 
codimY2 = di + d2. Then for — ^2 < M2, 

IM,M.^Y^^Y2) C iff + /f +-(^2), (5.8) 

2 ' 2 2 ' 2 
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while for M2 < —d2 

jM,,M2^Y„ Y2) C lf\{Y,) + C^-{Y2), (5.9) 



2 ' 2 



where we continue to use conormal notation. Here, the spaces on the left 

^2' 2^' 



hand side are conormal distributions with symbols of type 



We recall a version of the Cotlar-Knapp-Stein lemma [3l| that will be used 
in proving our Sobolev estimates. 

Lemma 11. Let j G V , T = "^j^^rTj with {Tjjjg^r bounded sequence of 
operators on , and {7(j)}jeZ'- be a multiparameter sequence of positive 
numbers such that A = J^jei.'- ^(^) ^ V 

\\t:t,\\ < {^{^ - j)r 

mT;\\ < (7(^-j))' 

then 

\\T\\ < A. 

Theorem 12. Let Mq be the operator that multiplies by q & where 

2 ' 2 

S has codim= d in X . For M = —d + a with < a < d, Mq maps to 

Remark 4. The assumption that a < d helps avoid multiplication by a 
distribution that has as strong of a singularity as the delta function. 

Proof: As Sobolev spaces are diffeomorphism invariant, we can apply an 
elliptic Fourier integral operator associated to a canonical transformation 
that turns q into 

j e''''-^'a{x-e')d9' (5.10) 

where a G 1 (l2|. 

2 ' 2 

Notice the kernel of Mq is equal to q{y)6{x — y). An interesting feature 
of this operator is that the kernel is also equal to q{x)6{x — y). Also note 
M* corresponds to multiplication by q. We take all these facts into account 
when composing the two operators. 
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Let {xj : j > 0} be a non-homogeneous dyadic partition of unity on M 
such that 

supp (xo)c{|t|<2} 

supp (x,) C {2^-1 < \t\ < 2^+1}, J > 1, and 
\xT\< Cm2~"'\ for all m > 1. 

In the language of Lemma [111 the terms Mg^ ^ will have Schwarz kernels 

^^^^-'^■^^'■\m)ciA^:d)dddi. (5.11) 



Set Xj,j',k,k' = Xj{\G\)x)j'i\T\)Xki\^\)Xk'i\v\)- Since the gj, ^'s are C°° functions, 
the compositions Mq_^ t^q , y well-defined and equal to M*, ^Mq., ^, as the 
multiplication operators are normal. The Schwarz kernels for these operators 
are 



Kjj,^k,k'ix,y) = j ^i[{x-zyii+z'-e-{y-zyr,~z'-r] ^ (5^^2) 

Xkil^DXk'ilvDO'jiz; e)aj>{z; T)dzdedTd^dr] 

where aj = Xj " « and aj/ = Xj' " o,- The numbers 7(j) will be of size as 
a result of our partition of unity. This kind of decomposition of an operator 
is commonly referred to as Littlewood-Paley decomposition; see jSlj. For 
brevity in our upcoming calculations, we refer to the phase in fl5.12p as $. 

Consider the operator = I — A^, where denotes the Laplacian in 
the 2;-variable. It follows that L* = L. Notice 



((e+(^,0)-r^-(r,0))) 



-e = e 



with the expression on the bottom being (Vx$)^^. It is the size of this 
Japanese bracket that will determine how we get our estimate in each region 
of integration. There are several cases we must consider: 

Case 1: Single elliptic variable. 

This region of integration corresponds to the part of the operator that is 
microlocally supported far away from the diagonal. Without loss of general- 
ity, because of symmetry between C, and r (likewise for 6 and r), we can let 
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the elliptic variable be ^. In this region, (V^$) ~ (^). Remembering that 
a spatial differentiation of symbols of type (|, ^) loses {O^y integration by 
parts Nz times with our differential operator Lz gives 



< J Xkm)xA\v\)xm)xA\r\) < e >"'^^< e >'^< e dzdOdrd^d^ 

< 2MU+j')+dU+j')+n{k+k)-N.k ^ (support of 2). 

(5.13) 

The above estimate is still rough as we have not utilized the integration in z 
and the later integrations in x or ?/ demanded by the application of Young's 
inequality; see 

If we let and be defined analogously to Lz, with the differentiations 
happening in the variables determined by the subscripts, then integration by 
parts and iV^ times in ^ and rj, respectively, gives us 

y e** X {2\x-z))-''^{Ljf^{2'''{y-z))-''^{L^f^ X 

(V.<l')-^^(Lf)^^[Xfc(|e|)Xfc'(l^l)%(^;^)%'(^;^)]t^-st^^t^^^^ec^^ 

Remembering (Lp^e(Lp^'' + (6,0) - r] - (t,0))-^- < {^ + {6,0) -r]- 
(r, 0))~^^ we get a better size estimate of 

where h = J (2^(x - z))^^^^ {2''' (y - z)y^^'^dz. Since i|x - y| < |x - 2] or 
\\x — y\ < \y — z\, it follows 

h<h = l{2'-\x~y))-'''H2''{y-z))-'''^dz + 

J {2\x - z))~^^i{2'''~\x - y)A^^"dz. (5.15) 

Using this inequality and completing the z integration in I2 followed by an 
integration in x or y, we see 

\\M,,',k,k'Ml,,,,A2 < 2^^(^-+^")+^(^-+^")-(2^^)^ (5.16) 
28 



Letting iV^ > M + (i, we get rapid decay of the norm. 

This kind of behavior is typical in regions far away from the diagonal of 
phase space. The following cases will show that such freedom for M is not 
allowed when phase variables are closer to the diagonal. 

Case 2 : > |r;|, |^| ^ |r|} or « |r;|, \e\ > |r|} 

Suppose we are in the first region of integration, (V^^) — v) ~ (0- 

The situation when ^ is elliptic to ^ or r is covered in Case 1. Otherwise, 
the gain obtained by (Vz^*) does not overwhelm the loss of (9)^ (r)^ . The 
reasoning for the second region of integration is identical. Hence, the upper 
bound for Kjj'^k,k'{x,y) becomes 

h X 2-^(^'+^'')+'^(-''+-''')+"('=+^') (5-17) 

Our only option left is to use the pseudodifferential operator part of Mq 
to get a gain in k (or k') after integration in z followed by an integration in 
a; or 1/ to get the additional gain in k' (or k) and finish the application of 
Young's inequality, just as in Case 1. This leaves us setting M < —d in order 
to sum along these indices. 

Because of the aforementioned symmetry in Case 1, the analysis in this 
region covers the situation when the elliptic variable in a pair is switched. 

Case 3 : {|^| ^ |^|, |^| |t|} 

The quantity (V^^) is bounded above by 1 as 

(V,$) ^{C + {9,0)-rj- (r, 0)) ^ {C - v) - h 

implying that the only gain we get comes from the pseudodifferential operator 
part of Mq, just like in Case 2. Following those arguments again, we must 
assume that M < —d. 

Case 4 : {|^| \r]\ ^ \9\ ^ \r\}. 

We are now on the diagonal and must estimate the operators Mjj whose 
kernels are of the form 

Since there are no spatial integrations of which to take advantage, we can 
only integrate by parts A^^ times in ^, giving us the estimate 
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Integrating in x or y gives an additional gain of 2 "-^ , forcing us again to have 
M < —d for summability. 

Case 5 : All other regions. 

It worth noting that outside of Case 1, all of the analysis done was dependent 
upon whether the gain obtained from the spatial gradient V overwhelmed 
the loss from the symbol. In each case, integration by parts in the phase 
variables of the pseudodifferential part of the kernel and the last spatial 
integration m x oi y always left us with 

\K,,,kA^,y)\ < h X 2*'^(^-+^")+'^(^-+^-')+"('=+'='), (5.18) 

before the z integration. Hence, all of the above analysis, including in Case 1, 
could have been done this way. The effort was made in Case 1 to show that 
microlocalization far away from the diagonal gives rapid decay, as expected. 

Therefore, by following the steps in Case 2 or after, any other region 
of integration requires us to have M < —d like in the previous cases. If 
M = —d + a, we see the norms grow at a rate of 2^°'. For a fixed e > 0, 
dividing out by 2^°''^^ gives the H~°'~^ norm of M„. , M* . ■ 

Corollary 13. As in Proposition [TOl with a < di, we have Mq : H^{X) — >■ 
H'-" for 

Ml + ^ < -di - da + a, for Ms > -^2 and Mi < -rfi - f + a 
Ml + ^ < -di - (ia + a, for Ms < -cis and Mi < -di + a 

Proof : Applying the parabolic cutoff to Mg, we get the sum Mg^ + Mg^ with 
qi and gs being of type (|, |). We apply our last proposition to obtain the 
necessary restriction on our orders. ■ 

Remark 5. Notice that for M2 sufficiently negative, Mi is allowed to be 
greater than —di. By estimates on the growth rate of nested conormal dis- 
tributions in Section this regime of orders allows q to blow up, i.e. be 
unbounded near Si. Moreover, using the fact that 

fl /^'^^'^HFi,F2) = /''nn), 
M2eM 

we have generalized an estimate used in f^J in the case of multiplication op- 
erators. This theorem will be stated in the next section. 
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6. Parametrix to the wave equation 



Parametrices for operators P{x, D) of real principal type, such as □ 



— A, were originally considered by Duistermaat and Hormander [5[ and 



were shown by Melrose and Uhlmann [17| to be operators whose Schwartz 
kernels are paired Lagrangians, whose wave fronts are contained inside the 
union of the diagonal relation and the flowout from the diagonal under the 
Hamiltonian vector field of the principal symbol of P{x, D). The fundamental 
solution D"^, which is an exact parametrix, lies in the space J^a'^a (A, Co), 
where 

Cn = {{x,t,uj;^,\^\,n),{x + Sj^^,s,u;^,\^\,n) : 

{x, t, w) G M" X M X ^ G M"\0, Q G T*{S''-^)} (6.1) 

is the canonical relation obtained by fiowing out from the light cone {|^| = 
|r|} along the vector field = ^ ■ ^ + 

The Schwartz kernel of D"^ can be written in the form 

J e'[{--y)<+{t-sM\+{t-s)p]^^^^ ^. p^^^^p (g 2) 

where the function appearing in the exponential is a multiphase parametriz- 
ing the associated canonical relations and a G S'~2'~2 (]R"+i x MP^^ x (M" x 
]R)\0). This generalized Fourier integral operator appeared in [51] but was 
not systematically treated until the work [17] and Guillemin and Uhlmann 
[ll| . Due to the presence of the diagonal relation, which is the wavefront set 
of standard pseudodifferential operators, operators of this type are given the 
name of "pseudodifferential operators with singular symbols"; an in-depth 
study on operators of this type is done in [9|. We will recall some map- 
ping properties of such operators, both on Lagrangian-type distributions and 
Sobolev spaces, from this article which will be of use when solving our inverse 
problem. 

The study of operators similar to from a microlocal perspective has 
been carried out in various geometrical settings and can more recently be 
found, for example, in the work of Baskin on scattering theory on de Sitter 



space and anti-de Sitter space. In [IJ], Joshi develops a symbolic calculus 



similar to that of Melrose and Uhlmann to construct complex powers of □ on 
Riemannian manifolds and develop parametrices for the resulting operators. 
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6.1. Action on spaces of distributions 

We now consider the mapping properties of a parametrix for a pseu- 
dodifferential operator of real principal type on various kinds of Lagrangian 
distributions. Of course, our intended application is for the d'Alembertian 
on M""^ X 5*""^, but because of the complexity of certain phase functions in 
our coordinate system, we prove the results in more generality in order to 
utilize the more computationally convenient normal forms. 

Let D) be an mth order pseudo differential operator with real homo- 
geneous principal symbol of classical type prn{x,^). If Cp = {(x;^), {y;rj) : 
(x;^) G char{P),(y,rj) G '^{x-,^} where '^{x-^) is a bicharacteristic of P{x,D), 
it follows from [l7| that the parametrix Q for P{x, D) lies in the class 
J3~™'~I(A, Cp). The following results are from [lo| . 

Proposition 14. Suppose Aq C T*{X)\0 is a conic Lagrangian intersecting 
char{P) transver sally and such that each ^characteristic of P{x,D) inter- 
sects Ao a finite number of times. Then, if T E P'\A, Cp), 

T:r(Ao)^r+^''(Ao,Ai), 

where Ai = Cp o Aq is the flowout from Aq on char{P). 

Proposition 15. Suppose Ai C T*{X)\0 is a conic Lagrangian which is 
characteristic for P, meaning that Ai C char{P). Then if T is as above, 
then 

T: r(Ai) ^r+P+^(A^) 

and therefore 

Q:r(Ai)^r+P-^(A,). 

Both of these propositions will be used to describe the scattering kernel when 
solving the direct problem. 

Because of the geometry that arises when dealing with an approximation 
to the scattering kernel, it is necessary to understand the action of an FIO as- 
sociated to Cp on paired Lagrangian distributions with wavefront containing 
characteristic points. The following theorem appears as Proposition 4.1 in 



111 in a more general form. We give an alternate proof in order to emphasize 



the multiphase interpretation of paired Lagrangian distributions. 

Theorem 16. Let u G /^■'^^'^^^(Ai, A2), with the Lagrangians intersecting in 
codimension d, and F G I^iT), where T is a homogeneous canonical relation 
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such that r o Ai = Ai and F o A2 = A2 and the compositions are transversal. 
If Ai and A2 are cleanly intersecting in codimension d as well, then 

e /*^i+P'*^2(Ai,A2). (6.3) 



Proof : By a result in [17|, we can conjugate by elliptic Fourier integral oper- 
ators associated with a canonical transformation and therefore assume that 
our domain is R", Ai = N*{{x' = 0}), A2 = N*{{x' = x" = 0}), where 
{x',x") G M™'"^"'. The oscillatory representation of u is now 

I e^ly'<'+y"<'%{y;e,ndede', (6.4) 

where b G S^^''^'^-^{W x (R""'^ x R'^)\0). For F, the Schwarz kernel takes 
the form 

/e-«.-%,.;«).. 

with (/)(x, y; 6) parametrizing our (transformed) canonical relation F. In order 
to show that Fu{x) is another paired Lagrangian, we break our region of 
integration into subregions and show the new phase is a multiphase and 
the product of symbols satisfies a symbol-valued symbol estimate in each 
subregion. The oscillatory representation of Fu{x) is 

j y. e)},i^y. ^")dedi'di"dy. (6.5) 

There are three regions to consider: 

Case 1 : If {|^'| ^ |^"|}, then u G I^^^^^^{K2)- The symbol h in (ES]) now 
satisfies a standard symbol estimate and we follow the proof of Hormander's 
result on the composition of Fourier integral operators whose canonical rela- 
tions intersect transversally. This requires us to define a new phase variable 
w = (d^P + 1^1^)2?/, 6,^)). It follows that in the region {|^| ^ \9\}, our phase 
in (16.51) parametrizes F o A2 and the product of a and h satisfies a standard 
symbol estimate. When {16*1 < |^'| ~ |^"|}, or vice versa, integration by parts 
in (16. 5p shows the order of the symbols can be decreased arbitrarily, making 
Fu a C°° function. Hence, Fu G /^i+^2+p(A2). 

Case 2 : If < then u G I^'+^^^{k2) once again and GG 

jj\/i+M2+p^yY^'j the above arguments. 
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Case 3 : Suppose {|^'| > |^"|}. In this region, we are microlocalized to the 
intersection Ai fl A2. First, if {y" 7^ 0}, then we can integrate by parts 
in (16 ■4p in C," arbitrarily many times to integrate out this variable and see 
that u e I^^'{Ai). Once again following Hormander's argument and setting 
to' = ((|^|2 + ICl^y^y^e,^')), we get that Fu G I^'^'+p{Ai). 

By our computations, it is natural to consider the new "big" variable in 
the paired Lagrangian distribution as u' = + \C,'\'^)^y,0,^')) and the 

new "small" variable as u" = Let us recall that on Ai, 

dy(f){x,y;e) + {C,0) = 
d0(f){x,y;9) = 

y' = o, 

and on A2, 

dy<P{x,y;9) + iC,n = Q 
dg4>{x,y; 9) = 

y' = y" = 0. 

Hence, in a conic neighborhood of the intersection, we are near the closed set 

(dy(l){x,y-9) + {^',0) = 

lde<f){x,y;9) = (6.6) 
[y' = y" = 0. 

As the phase function {/)(x, y; 9) is homogenous of degree 1 in the 9 variable, 
we can assume \dy(f)\ ~ 16*1. It follows that within a conic neighborhood of 
(I6.6p . 1^1 ~ |^'|. The symbol in (I6.5p . which we will refer to as c{x;u',u"), 
will now lie inside of S^^+P'^\ ■ 

Remark 6. A corollary to the version of this theorem that appears in fl2] 
says ifV is a homogeneous canonical relation, satisfying the same transver- 
sality conditions with respect to the codimension d intersecting pair of La- 
grangians (Ai,A2), then (F o Ai,F o A2) is a codimension d cleanly inter- 
secting pair once again. We could have used this theorem to do away with 
the clean intersection calculations in Section\^ However, the computations 
themselves give us a deeper understanding of the present geometry, reappear- 
ing in calculations that show the original surfaces Si and S2 are one-to-one 
with our flowouts and reflected Lagrangians, which is the crux of the inverse 
scattering problem. 
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Let us now analyze the application of D"^ to a nested conormal distribu- 
tion that will appear in a later section. Using the oscillatory representation 
of q5 in I'^^^'^'^^Si^, S2+) discussed in Section [5l 

□-1^^ = j ^i[(x~y)<Ht-sM\Ht-s)p]Hs-y-'^)r+Myy9'+h2{yye''] 

xa(x, y- p, i)b{y- 6', 9") dpd^dTde'de"dyds. (6.7) 

where b is localized to the region {|r| ^ \6'\ > \6"\} in order capture the 
intersection Ai+ fl A2+. 

To microlocalize to Cn H A, we localize a to the region {|^| > If we 
do not focus on this part of phase space, then we are microlocalized to A 
away from Cn and are therefore applying a pseudo differential operator; the 
wavefront set of D~^q6(t — x ■ uj) does not move if D^^ acts like a pseudodif- 
ferential operator. Wavefront set calculus tells us 

WF{U-\5{t -x-u))c. Ai+ U A2+ U A^+ U A^+. (6.8) 

The calculations in Section H] show that these Lagrangians all share a common 
submanifold of intersection, namely Si U S2. By interpreting the phase in 
(16. 7p as a generalized multiphase, integration by parts and application of 
the standard theory of Fourier integral operators in (16. 7p shows that these 
singularities do appear. A theory of intersecting quadruples of Lagrangians in 
a certain configuration does appear in [17] but does not apply in this setting 
because certain conditions involving the arrangement of the Lagrangians are 
not met. 

If we further localize our symbols to {|p| ~ |r|}, then letting our "big" 
variable be (d^P + + )2 (y, s), ^, r, ^0, 

our "medium" variable be p, 
and our "small" variable be 9" , we get a distribution whose wavefront set is 
now inside of Af^ U Ai+ U A2+. We will elaborate further on this triple in 
Section [71 

The final microlocalization will be to further localize our symbols to the 
region {|^| ^ |r|} in phase space and to {t 7^ s} in space. Integrate by 
parts arbitrarily many times in p to obtain a new symbol c(x, y\ ^, r, 6'', 9") = 
b(y;9',9") x J e^^^^~^^''^a{x,y;$,,p) dp satisfying a symbol-valued symbol esti- 
mate 

d2,ydS.dl^^,,cix,y;^,T,9',9") < (e, r, 0', ^^'V^-^-H ^ (6.9) 
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We now let our "big" variable be ((|^P+|r|2 + |0'|2) 2 (y, s), ^, r, ^0 and "small" 
variable be 9" . This stratification puts the distribution into the regime of 
a paired Lagrangian distribution associated to the codimension ^2 cleanly 
intersecting pair of (A^_,_, Ag.,.). Section [8] shows that this portion of 
contains all the information needed to solve the inverse problem. 

6.2. Sobolev estimates 

We equip T*(M"')\0 with the canonical symplectic two-form u = TidC,i A dxi. 
Let E C T*(R") is a smooth, codimension k conic submanifold with 1 < k < 
n that is also involutive, meaning that the ideal of smooth functions that 
vanish on E is closed under the Poisson bracket. Theorems in symplectic 
geometry from Section 3 of 0] show that the flowout of S, As G T*(M"')\0 x 
T*(M"')\0, is a canonical relation given by 

Ce = {{{x;0,{y,v)) G S X S : G E^,,^)}, (6.10) 

where is the bicharacteristic leaf of S containing (x;^). It is straight- 

forward to see Cs and A intersect in codimension k; refer to Section 4 of [^. 
The following is a theorem from j^. 

Theorem 17. Let A G P^^W x W; A, C^). Then 
continuously for a// s G M if 

k 

sup{p+ + < -So. 

Consider the flowout relation for D"^ in (16. ip . The parametrization shows 
that the cotangent variables to {x, t) sit in a hypersurface in the fibers of 
our phase space, namely the light cone. Because this submanifold has codi- 
mension 1, it follows by a rank calculation that A; = 1 in the above theorem 
for A = n^^. Hence for — sq > max(— 1, —2) = —1, O'^ is smoothing of 
order at least 1. In fact, away from the characteristic variety, D"^ acts as a 
pseudodifferential operator of order —2 and increases Sobolev regularity by 
2. 

Remark 7. It turns out that multiplication by standard conormal distribu- 
tions fall under the hypothesis of this theorem. 
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7. Born series 



We analyze the formal series 

oo 

j2{-ma-'M,mt-x-uj)), (7.1) 

i 

which in the physics literature is called the Born series , in order to solve the 
direct problem. Set Ui := {—iy{D~^Mg){6(t — x ■ u)) where Uq = 6(t — x ■ u) 
so that the series on the right of (17. ip is formally telescoping when □ + g is 
applied. The remaining parts of this paper will show the first two terms, 

UQ + ui=5{t-x-uj)- (□"^Mg)((5(t - X ■ w)), 

commonly referred to as the (first order) Born approximation, will be enough 
to solve the inverse problem. Analysis of the higher terms is substantially 
more intricate, involving a more complete theory of distributions associated 
to higher order systems of Lagrangians. At this point, we do not know if an 
in-depth analysis is even possible without strict assumptions on the geometry 
of Si and 5*2. 

7.1. Computing the series 

For q G /^^'^^^(S*!, 5*2), we apply a microlocal partition of unity Yl^=i Xii^j ^? ^) 
that gives 

xiq e /*^H5i), 

X2q e I^''+^'%S2), and 

Xsq e /*^-^'^^(^i,^2), (7.2) 

where the last term is micro locally supported near the intersection Ai H A2. 
This partition is introduced to facilitate the analysis of the singularities of 

Ui. 

By Lemmas [H] and [9] and remembering that supp 6 G {t — x ■ u = 0}, 

Xiq-6it-x-cu) e /°'^^H5+,5i+), 
X2q-S{t-x-u) e /°'^^i+^^^(5+,52+), and 

X3q-S{t-x-u) e J°'^i'*^H5+,Si+,52+) + /*^-"''^(Si+,^2+). (7.3) 
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Keeping in mind the "good" part of U T? does not meet A_|_ fl Ai_|_ and 
A+ n A2+ , or rather A+ n Ai+ C S^\Si and A+ n A2+ C E^yEs, we can mi- 
crolocahze away from a conic neighborhood O of (S^\Si)U(S^\S2), changing 
(Q to 

Xiq ■5{t-x-cu)e /^-?(A+\(0 n A+)) + /^^i+^-t(Ai+\(0 n Ai+)) 

1 n 

X2q ■6{t-x-u)e r^'^{A+\{0 n A+)) 

+ I'^^^^'^^'^-HA^AiO n A2+)) + H-\ 

1 n 

Xsq ■5{t-x-uj)e l2-^{A+\{0 n A+)) 

+ jM.+M.+i^l±^-t,-M.-^ (^^^^(0 ^ Ai+), A2+\((^ n A2+)) 

+ H^^\ (7.4) 

where we have switched from conormal notation to Lagrangian notation and 
H'^ are elements w G if*"(]R" x M x S"""^) with WF{w) C O. Note the sums 
of spaces are not direct sums as there is some overlap in the singularities 
from one term to the next. 

We are now ready to describe the second term, ui, of the Born approx- 
imation. Let L := Co o O, where Cn is the flowout relation for D"^. The 
set L will be invariant under the Hamiltonian flow of H^. Applying D"^ to 
both sides of f l7.4p . using Propositions [T4l and ITSl and Theorem [T6] gives 

u, = n-'q6 G /^-i(A+\L+) + J*^^+^-t-f'-i(Ai^\Li+,A?+\L?+) 

+ /^■^^^^^^^^-5-i'-nA2A^2+,A^,^\L^^) 

+ /(A?+\L?+, A^+\L2^, AiALi+, A2+\L2+) + i/f +^ (7.5) 

where the last term I{Al^\Ll^, A^+XL^^, Ai+\Li+, A2+\L2+) in (1731) is a 
new class of distributions such that WF{u) C A5^_^UA^_^UAi+UA2+, H^j^'^^ are 
Sobolev elements of order sq + I with wavefront set inside of L, and the varied 
labeling of L denotes intersection of L with the juxtaposed Lagrangian. Note 
that we have used Theorem [T7] on the Sobolev space term. The quadruple 
portion of Ui was analyzed earlier in Section [61 

In Section [6l we prove, microlocally away from A"_^ fl A2+ fl Ai+ fl A2+, 
there exists a nontrivial singularity on A"^ fl Ai+ fl A2+. We label this class 
similarly as /(A"^, Ai+, A2+), the space of distributions whose wavefront set 
is contained inside the union of the three listed Lagrangians. Notice we do 
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not attach orders to this newly defined class in spite of it being similarly 
defined to those of the nested triple conormal configuration. This is because 
for codimensions of 5*1 greater than 1, there exists a conic singularity on 
along 5*1+. Hence, because of the degeneracy of the underlying submanifolds, 
this new triple does not fall into our previously defined class. However, for 
di = 1, the distribution is a nested triple conormal of order Mi, M2). 
The space 

also appear as a portion in the description of the quadruple term from (17. 5p . 
After all these observations, it follows that 

+ /ioc(A?+\L, A^+\L, Ai+\Li+, A2+\L2+) 
+ /ioc(A?+\L?+,Ai+\Li+,A2+\L2+) 

where the subscript "loc" denotes localization near the respective triple and 
quadruple intersections. As the t variable is bounded on Ai+ and A2+ because 
the support of ui is inside {t — x ■ u = 0}, for t » 0, 

ui e /^-^(AAi^+) + /*''^+'^-^-^(A?A^?+) 
+ /^'^^+^'^^+^^-t-i(A^+\L^+) 

^ ^M,+M2+i^i±|^-t-i-^2-f (A«^\L^^, A^2+\^2+) + (7.8) 
7.2. Subsequent terms and singularities 

We show in this section that the remaining part of the Born series, i.e. u = 
U2 + + has WF which grows larger with each iteration of □~^Mg. 
As a consequence, we will use the mapping properties of D"^ and Mq from 
previous sections to place the u into a suitable Sobolev of higher regularity 
than the previous two terms. This final step will be taken in Section [HI 

It was shown in Section H] that A"^ and A2+ are conormal bundles of 
the hypersurfaces S^^ and 82^, respectively. A well-known fact about the 
multiplication of distributions u,v & V{X), where X is an open set in M"', 
states if WF{u) n WF'{v) = then 

WF{uv) C WF{u) + WF{v). 
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wheTeWF{u)+WF{v) = {ix,Cu+Q : e WF{u),{x,C,) G WF{v)} C 

T*(M" X R X S"-^)\0; see |12|. 

If we microlocalize u G /^i+^2+ ^ 2^ 5"I'^I(A2+, A2_|.) away from A2+, 

we get u G IMi+M2+^^ t-2(A^^). In fact, u is conormal to S'g.,. away 

from A2+ n I\.2j^. It is the interaction of this part of ui with xi ^ /*''^^(>S'i) 
that creates the new WF. 

Using Lemma [H all we must show is that there exist submanifolds 5*1 
and 5*2 such that Si iti 6*2+ and this intersection is different from all other 
submanifolds of the same codimension. 



Proposition 18. Let Si = {x„ = 0} and S2 = {xn = Xn-i = 0} be subsets 
jjn+i ^ Then Si ftl S!^^ is a new submanifold of codimension 2 in 

W^^^ X E"-^^ that is different from those that ui is conormal. 

Proof: Considering the general parametrization of A.2+ given in Section 5, we 
apply the Implicit Function Theorem to solve for the yn-2 coordinate. More 
explicitly, 

:i/„_2 = /(l/i,...,l/„-3),(^i,^2) eM2\0, 

uj G §"-\r G M} 
From the parametrization of Si, which is 

{((l/",Z/n-i,0),s,u;) : (i/",y„_i) G M"-\ s G M,u; G S"-^}, 

it follows 53 = Si n 5*2+ has the same parametrization as 5*2+ but with 
~ 2{eiLd^\'+e2Uj ) ' 'k ' "^"^^ brevity in the following argument, set r = 



9'f+e? 



2 



2{dlLUn-l+d2UJn) 

Computing tangent vectors, we see {{■^}7=i^'§;^{'^}7=i} spans 



a 2n — 1 dimensional set in Tx{W^^^ x ^) for x G 5*3. Take ^ on 5*2+ and 

observe its nth spatial coordinate, which is ~ 'Ml'^^'n)- This quantity is 
not identically equal to because 0-sections are deleted from the total space 
and the ability to work away from (g|^r)w„ = 1, a lower-dimensional set that 
is independent of r. Hence, at a; G ^3, T^{Si) © = T,.(M"+i x S"-^) 

and the intersection is transverse. 
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We now show 5*3 is not one of the previously encountered submanifolds. 
It is clear 5*2 has its coordinate identically equal to 0. However, -^{Oi — 
TUJn^i) 7^ 0, almost everywhere. This shows S2P\S^ C S^. In 51+, the relation 
t — a; ■ a; = is satisfied. When this relation is applied to the parametrization 
of 5*3, it follows that —r{6iUJn-i + 020Jn) = 0. This is another equation that 
holds on a lower dimensional set, proving 5*1+ H 6*3 C 5*3. ■ 

It follows from Lemma [H]that q{x) ■ ■ 5{t — x ■ u) has different 

wave front set than does uq + ui. This gives the evidence, after application 
of D"^ to the above expression, that the new fiowouts generated are different 
from A"_,_ and A2_,_, as N*{S3) has characteristic points because 5*3 C S'f^. 
Hence, the wave front set of the Born series does not stabilize. We will show 
in the next section that the approximate solution uq + Ui will be enough to 
solve the inverse problem by simply including the remaining terms of the 
Born series in a Sobolev space. 



8. Solution to the inverse problem 

We now solve the inverse scattering problem by combining all of the results 
from previous Sections. First, we define the scattering kernel, which describes 
the far-field pattern of a solution u to the direct problem, using the theory 
of Lax and Phillips [l6|. We then show that the restriction of the scattering 
kernel, known as the backscattering, also determines the singularities of q. In 
fact, we prove a corollary that shows for well-behaved subsets of scattering 
data D, the corresponding restriction of the scattering kernel continues to 
determine the singularities of q. 

8.1. Lax-Phillips theory and the scattering kernel 

The Lax-Phillips Radon transform 16| is the map Rlp '■ P'(C^) — t- V'iG) 



and is defined by 

Rlp{v^,vi) = Cn\D;^\{D,R{vo) - R{vi)) (8.1) 



where C„ is a constant depending on n, \Ds^ | is a pseudodifferential oper- 
ator acting on the s variable of order which overlaps with a differential 
operator for n odd, and R is the Radon transform defined in Section HI To 

n — 3 

define the scattering kernel, we set w = ^ \{u{x, t, u) — 5{t — x ■ u)) and 
get the relation 

aq{t - s, 0, co) = S{t - s) (g) 5{(f) -co) + Rlp{w, Dtw), for t » 0. (8.2) 
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Focusing our attention on the scattered part of u which is expressed as m — 5, 
it follows that 

{aq-6®6){t- s,(p,uj) = F{u-6) (8.3) 



-n-l 



where F G I^^{Cr) and Cr is the canonical relation associated to the 
Radon transform; see Section HI Since F is elliptic on the first components 
of CRoh.1j_ and CroK^j^, these are the only sets that stay after the application 
of Cp, also defined in Section HI The reason for introducing p* is to deal with 
the translation invariance of fl8.3p for t » 0. Summing all this up, we get 

Definition 16. The scattering kernel associated with q is defined by 

{ag-6^6){s,(f),u) = p*F{u-6), (8.4) 
where u is a solution to the continuation problem. 

Unlike in the case for scattering theory for the wave equation when q G 



C^(M") where there is a wave group [18[, we use the Born series to provide 



information about the solution u. It is natural to consider the approximate 
scattering kernel p*F{uo + Ui — S)= p*F{ui) and investigate the error of the 
difference between this and the true scattering kernel. The next section will 
show this error term is of lower order, giving us the ability to work with the 
approximate scattering kernel instead of the exact one. 

8.2. Comparison of scattering kernels 

Our principal objective is to make a comparison between aq and p* F{uq+ui — 
(5), the exact scattering kernel and the approximate scattering kernel. All of 
the analysis from our previous sections and further computations below will 
show that this difference of scattering kernels is negligible in an appropriate 
sense. 

For now, set u := u — [uq + Mi). For t » 0, 

(□ + q)u = — Duq — q ■ Uq — Dui — q ■ Ui 

= -q-ui (8.5) 

We recall that for the same range of t, 

/*'^^+^-t-f(A^^\L?^) 
-i(A^+\L^J 

-i-M.-^(A«^y^a^^A^2+\^24-) + (8-6) 



Ml G 


T — 
I 2 




+ 




-M2+ "1+^2 + 1 n 


+ 


jMi4 
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Moreover, 

p*F{u-5)- p*F{uo + Ui-5) = p*F{u- (uo + ui)) 

= P*Fiu), (8.7) 

showing why placing p*F{u) into a space of lower order will allow us to ex- 
tract information from the true scattering kernel by reading the approximate 
scattering kernel. It is necessary to prove this error term does not interfere 
by computing it's Sobolev regularity. 
We define 

oo 

(/ + n-'M^r' := ^(-l)^(n-iM,)^ (8.8) 

j=0 

as this makes the Born series consistent with a Neumann series; also, we 
only need to know the Born series asymptotically. In addition, □~^Mg must 
raise regularity on Sobolev spaces in order to have the subsequent terms of 
the Born series be smoother. This happens when the orders of q satisfy 
Ml + ^ < -di- d2 + l for the specific ranges of Mi and M2 discussed 
in Theorem [121 Recall that this allows q to blowup on parts of Si or S2 
depending on the specific orders. Applying D^^ to both sides of (18. 5p gives 
the new relation 

{I + D-^Mg)u = U2:=n-^Mg{ui). (8.9) 

Because of the iterated regularity characterizations of the spaces that appear 
in (18.61) . we can place Ui into a Sobolev space of some order s. Another 
application of ^^Sf^\-iy {D-^ M^Y to both sides of ([83]) leads to 

N-l 

{-i)^-\n-'Mgfu = J2i-^yio-'Mgy{u2), (8.10) 

j=0 

with n-^Mq raising Sobolev regularity by some amount in each term. It 
follows that the right hand side of (I8.10p stays in the Sobolev space H'^ by 
convergence of the Neumann series. Hence u G H^. Set L := CpO Cr o L. 
Keeping in mind the mapping properties of p* and F on Sobolev spaces 
[7I.I12I and spaces of Lagrangian distributions from the transverse intersection 
calculus discussed in Section [21 (18. 6p implies 
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Theorem 19. The difference of the scattering kernels, a{s,(f),uj) — p*F{uq + 
Ui — 6), is in and 

(8.11) 



G 




f) 




+ 




(A? 




+ 




3,3 
2 + 




+ 




1,3 
2^ 




+ 


„ 1 3-2n 







Hence, the principal symbol of a is same as that of p* F{uq+ui—5) = p*F{ui). 
8.3. Determination of Si, S2, and p{q) 

A majority of this final section is essentiaUy an adaptation of the conclusion 
of 10| to the present context, with some minor modifications. In order to 
even retrieve the singularities of q from the restriction of a{s, (p, u) to various 
sets of scattering data, i.e. submanifolds of M x 5*""^ x S*"^^, we must be 
able to do this without any restriction, i.e. just from using the full scattering 
kernel. By Theorem [19] of this section, it is enough to show Si and S2 are 
determined by the refiected Lagrangians. We will compare two refiected 
Lagrangians: 

{y,v) e N*{Si),cue S''-\r eR,a e M\0} (8.12) 

and 

Ai+ = {{-y{- — = — - + uj),-^ — - — '-,u-a,-m%_,^Ay),^Ciy)) ■ 

a a ^ 

(y,z/) G N*{Si),u G S"-\f G R,a = with z/ • d; ^ 0](8.13) 

2u ■ CO 

First noticing that co = u and a = a, we set the coordinates equal to get 
u = u. The comparison of the Q coordinates tells ns y = y + c{y, uj)ijJ. Now, 
set the s coordinates equal to get the relation 

a a 
^-y-iy = -y-iy^y-iy={y + c{y, u)u) ■ v 

^ c{y,uj)u ■ u = (8.14) 
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Remembering ■ z/ 7^ because of the tangential rays condition, we get 
Ci{y,uj) = 0. Therefore y = y. Hence AJ^ determines Si. The same holds 
true for A2+ and 82- 

Let B = {{s,(f),uj) eRx S"""^ x S""^ : = -u} be the backscattering 
surface. The map 

jB : K X 5""^ ^ 1 (8.15) 

defined by jM{s,uj) = {s,ijJ, —u) induces the pullback 

: P;(M X S"-^ X S"-^) ^ X (8.16) 

with the domain being distributions whose wavefront set is disjoint from 
the normals of jb- This pullback is another Fourier integral operator, e 
/^(Cb) with 

s = s' ,uj = uj' = —(f), (r, Q) = (r', Q') such that 
(r',$,fi')^Ar(;_,.)W} (8.17) 

We set Lb := Cb o L and view this as our new "bad" set. As A+ C N*{{s = 
0, = uj}) it follows that Cb o A+ = 0. It is easy to verify the compositions 
with our reflected Lagrangians are again transverse, yielding 

A^ = {(-2y.^,^-\u\,\u\zUy):{y,u)eN*{Si)} (8.18) 

= {(-2i/-^,|^;|H,kK^(y):(y,//)GiV*(S2)}. (8.19) 

A calculation similar to those in Section H] shows Aj^flAj^ is again another 
codimension ^2 submanifold of both Aj^ and A^. The polar coordinate argu- 
ment in the final part of Section H] proves this intersection is clean. We apply 
Theorem [16] from Section [6] to get the first part of our main result. 

Theorem 20. The backscattering, the full scattering kernel restricted to the 
backscattering data, is 

«|B = jB(a) e 

^-i+M.+M.+^-t-l+^.-M-^(AS\(LB n as),a^\(Lb n A^)) 
+ /^+^'^^+^-t-i+^(AS\(LBnAS)) 
+ /^+*^^+*'^^-^^-t-f+^(A^\(LBnA^)) 
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and determines the suhmanifolds Si and S2 as well as fi{q), for Sq and s as 
in Theorem [7P1 



Proof : Consider the parametrization of in f lS.lSp . Then ^ — ^su = 
'i*j^{y) + (y- ]^)(]^) = y- This means that Si is determined by A^. The same 
holds for 5*2 from the parametrization of Aj^. 

Theorem [19] tell us the principal symbol fi{a\s,) of the exact scattering 
kernel is equal to fi{p*Fn~^{q ■ 6)). The ellipticity of F,p*, and implies 
/i(ajB) is a non-zero factor times fi(D~^{q6)). By results on the symbol 
calculus of FIOs in [l2|, we are able to divide out by these multiples and 
focus on fi(D~^{q6)). 

Microlocally near Si U S2, acts as an elliptic FIO associated to Cn 
described in (16. ip . Recall Theorem [16] and the discussion after Remark [6] 
which sets up the oscillatory representation of the paired Lagrangian in f l8.6p . 
We see fi{D~^{q6)), in the coordinates given by f l4.15p or fl4.16p . is another 
non-zero elliptic factor times cr(g5) on A"_,_ U Ag^, by symbol calculus results 



As /i(n~^(g5))(x — r(z/ — au), x-uj + r,uj; u — au, a, — crz* (a;)) is a function 
whose variables parametrize the bicharactersitics that foliate A"^_ and A2^, 
we can flow our symbol back in the r to Si U S2, obtaining fi{n^^q6) = 
fi(D^^){x, X ■ co,u; u — au, a, — cri* (x)) x fi{q){x, u) x cr(5)(x, x ■ uj^uj] u, t, Q) 
near Si U S2. Note that fj,{6){x, t, u; u, t,Q) = 1 and fi{D~^){x, x ■ u — 
aoj,a, —ai^{x)) 7^ in this region. Therefore, we can divide out by these 
elliptic factors, leaving us with fi{q){x, u). As Si U S2 is a dense open subset 
of S^ U S^, it is possible to recover yU,(g) on all of Ai U A2 as is function of 
only X and z/, with no restrictions on either variable besides lying in Si U S2. 
■ 

It is not hard to generalize this result. In fact, the operator that 
restricts to a submanifold D is another elliptic Fourier integral operator that 
satisfies the necessary transversality conditions. If we have the Lagrangians 
associated to our new scattering data aijo cleanly intersecting in codimension 
(i2. Theorem [20] still holds. Since this last step follows immediately from the 
general form of Theorem [16] that appears in [ll| , we have 



in [nj. 





for all s G M, 
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3. if ^(s, z/) • z/ 7^ for all v in the image of the Gauss map of Si and S2, 
for every s G R, 

then a\o determines Si and S2 as well as fi{q). 

Remark 8. Conditions 1-3 are sufficient to carry out a calculation similar 
to that in the proof of Theorem ESI describing the steps to reconstruct the 
submanifolds Si and 5*2. 



9. Appendix: Multiphases 



The following is a proposition from [20 



Proposition 22. Let Ao,Ai be Lagrangian submanifolds ofT*{X)\0, Xq G 
Aq, Ai andpi be a homogeneous function of degree 1 such that pi{Xo) = and 
the Hamiltonian vector field Hp^ associated to pi is not tangent to Aq. //Ai is 
the fiowout from Ai fl {pi = 0} by Hp^, then there is a multiphase function ip 
that parametrizes the pair (Aq, Ai) which can be chosen such that ^{x, s, 6) = 
Pi{x, dx(p) and (p{x, 0, s) = ipo with ipQ a phase function parametrizing Aq. 

Remark 9. Proposition [H can be generalized to the situation where the 
Lagrangians cleanly intersect in a higher codimension. 

We will use this proposition to determine a multiphase functions for the 
pairs (Ai_|_, A"_,_) and (A2+,A2+) that appear in Section [71 Let us compute 
the phase for first pair (Ai_|_, A"^). By the proposition of Mendoza, obtaining 
the desired multiphase involves solving the following initial value problem, 

!f{x,t,uj;e,a,0) = (fo 

where t) = ^ — r defines the characteristic variety of □ in T*(R" x R x 
S*"""*^) and v5o(a;, t, uj\ 6, a) = H{x) ■ 9 + (t — x ■ ijj)a. We are subsequently led 
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to the system of ordinary differential equations 



j„ — _9 -L 



< 



dr 
dt 
dr 

dk 
dr 
dr 
dr 
ds 
dr 



2^, 

2 



for i = 1, ...,n 
for i = 1, n 



0, 


1 

s(0) = 

Solving for the characteristics, after replacing r for s, we come to 



(9.2) 



Xi{s) 

tis) 



-2—)s + x°, for i = 1, ...,n 

T 



(9.3) 
(9.4) 



Now solving for x° and t° in terms of s and plugging the resulting equations 
into ipo yields 

(f = Hi[x H s) - 9 + (t+ ( h l)s - (x s ■ uj))a (9.5) 

as our desired multiphase. Substituting H2, the defining functions for S2, for 
f^i in (19. ip and (19.51) gives us a multiphase for the pair (A2+, Ag.,.). 
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